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STRONG SCARRING OF LOGARITHMIC QUASIMODES 


SURESH ESWARATHASAN AND STEPHANE NONNENMACHER 


Abstract. We consider a semiclassical (pseudo)differential operator on a compact surface 
(M, g), such that the Hamiltonian flow generated by its principal symbol admits a hyperbolic 
periodic orbit 7 at some energy Eq. For any e > 0, we then explicitly construct families of 
quasimodes of this operator, satisfying an energy width of order g | lo g ^| in the semiclassical 
limit, but which still exhibit a “strong scar” on the orbit 7 , i.e. that these states have a 
positive weight in any microlocal neighbourhood of 7 . We pay attention to optimizing the 
constants involved in the estimates. This result generalizes a recent result of Brooks [Brl3] 
in the case of hyperbolic surfaces. Our construction, inspired by the works of Vergini et al. 
in the physics literature, relies on controlling the propagation of Gaussian wavepackets up 
to the Ehrenfest time. 


1. Introduction 

The main objective of Quantum Chaos is to understand the behavior of wave or quan¬ 
tum propagation when the underlying classical Hamiltonian dynamics is chaotic, or at least 
presents various forms of hyperbolicity. Semiclassical analysis shows that the wave mechani¬ 
cal system should “converge towards” its classical counterpart, in the high-frequency regime; 
the latter is equivalent with the semiclassical regime h —)■ 0, where h is an effective Planck’s 
constant. However, when the classical dynamics presents some hyperbolicity, this semiclas¬ 
sical correspondence breaks down when the time of evolution exceeds the so-called Ehrenfest 
time T e ~ c 0 | log /z|, where c 0 > 0 depends on the classical dynamics. As a result, a precise 
description of the stationary states (eigenstates) of the quantum system remains a challenge, 
since a precise description usually requires one to understand the quantum evolution up to 
much longer times. 

Although our results apply to more general situations, we will in this introduction focus 
on the case of the Laplace-Bcltrami A g on a compact Riemannian manifold ( M,g ). To 
describe the high frequency dynamics generated by this operator, we will use a semiclassical 
formalism, that is define the semiclassically rescaled Laplacian 

P(h) = -n 2 A g , 

and consider its spectrum in the vicinity of a fixed value E = 1, in the semiclassical regime 
h <C 1. The operator P(h) generates the dynamics of a free quantum particle on M, embodied 
in the Schrodinger equation 

= p {h)<p{t), (i-i) 

where ip(t) € L 2 (M,dg ) is the wavefunction of the particle at time t. In the semiclassical 
limit, this dynamics can be related with the free motion of a classical particle on M, namely 
the geodesic flow. This flow, which lives on the phase space T*M 3 (x, £), is generated by 
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the Hamiltonian function p(x,£) = ||£|| 2 , so we will denote it by By scaling, we may 
restrict this flow to the unit cotangent bundle S*M = p _ 1 (l). 

When M is compact without boundary and the metric g has everywhere negative sectional 
curvature, the geodesic flow is of Anosov type: every trajectory is hyperbolic. Ergodic theory 
shows that this flow is ergodic and mixing with respect to the Liouville measure on S*M. 
Anosov flows have been thoroughly studied, and represent the strongest form of chaos. One 
objective of Quantum Chaos is to study the localization properties of the eigenstates of P{h) 
at energies E h ~ 1, whilst taking into account the chaotic properties of the geodesic flow. 
The semiclassical regime h —> 0 amounts to studying the eigenmodes <pj of the Laplacian of 
frequencies Xj ~ k 1 —* oo. 

The Quantum Ergodicity Theorem [Sch74, CdV85, Zel87] shows that if the geodesic flow is 
ergodic, then almost all the high-frequency eigenstates p 3 are asymptotically equidistributed 
over M. Namely, there exists a subsequence dcNof density one, such that the probability 
measures \pj(x)\ 2 dx weakly converge to the normalized Lebesgue measure when S 3 j —$■ oo. 
An analogue of this theorem is established in the boundary setting for Dirichlet eigenstates 
by Gerard-Leichtnam [GL93] for domains in R n with W 2, °° regularity and Zelditch-Zworski 
[ZZ96] for compact manifolds with piecewise smooth boundaries. 

An open question concerns the existence of exceptional eigenstates, which would localize 
in a nonuniform way on M. Rudnick and Sarnak conjectured that for (M,g) of negative 
curvature, such exceptional eigenstates do not exist [RS94], a property called Quantum 
Unique Ergodicity. So far this conjecture was proved only in the case of certain hyperbolic 
surfaces with arithmetic properties [Lin06, BL14], Numerical computations of eigenstates of 
ergodic Euclidean billiards [H84] have shown that certain high frequency eigenstates present 
an enhanced intensity along short closed geodesics. In the case these orbits are hyperbolic, a 
quantitative characterization of these enhancements — baptized as scars of the geodesic 7 on 
the corresponding eigenstates — has remained elusive, in spite of intensive investigations in 
the physics literature (see e.g. the review [Kap99] and references therein), or a more recent 
numerical study by Barnett [Bar06]). In particular, it is unclear whether these enhancements 
can take the form of singular components w 7 5 7 , w 7 > 0 , in the weak limits of the measures 
\{pj(x)\ 2 dx. 

To connect the localization of eigenstates with the classical dynamics, which takes place in 
the phase space T*M, it is convenient to lift the localization properties to phase space, that 
is characterize the semiclassical microlocalization of the eigenstates. More generally, for any 
semiclassical sequence of L 2 -normalized states ((pn)n->- o> we may characterize the asymptotic 
phase space localization of these states through the semiclassical measures associated with 
this sequence. Let us recall how these measures are constructed. The sequence of normalized 
states (iph)h -^0 is said to converge towards the semiclassical measure p sc if, for any observable 
a G C™(T*M), one as 

(<Ph, Op s (a)^) / a dpi sc . 

Jt*m 

Here Op h (a) is an operator on L 2 (M), obtained through a semiclassical quantization of 
the observable a(x,£) (see section 2). The limit measure is independent of the choice of 
quantization. For any sequence ((pn) one can always extract a subsequence (ph,) converging 
to some semiclassical measure p sc . The latter is then said to be a semiclassical measure 
associated with the sequence (pn)- Semiclassical measures provide a notion of phase space 
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localization, or microlocalization : a sequence (<£>&) will be said to be microlocalized in a set 
K c T*M if any associated measure is supported inside K. 

Coming back to the scarring phenomenon, a sequence of eigenstates (tpn) of P(^), of 
energies E^m 1, is said to exhibit a strong scar on a closed orbit 7 G S* X if any associated 
semiclassical measure contains a component iu 7 <5 7 with some positive weight w 7 > 0 (here <5 7 
is the normalized, flow invariant Dirac measure on the orbit 7). Hassell has shown [HasslO] 
that such strong scars may exist in the case the orbit 7 is marginally stable (and belongs to 
a family of “bouncing ball orbits”). In negative curvature, the existence of such exceptional 
eigenstates would contradict the QUE conjecture. In [A 11 O 6 , AN07] it was shown that the 
eigenstates cannot fully localize along 7: for any sequence of eigenstates, the weight u> 7 
is necessarily smaller than unity, in particular in the case of constant negative curvature 
one must have w 7 < 1 / 2 . So far exceptional eigenstates have been exhibited only for a 
toy model of Anosov system, namely the hyperbolic automorphisms of the 2-dimensional 
torus, casually called “Arnold’s cat map” [FND03]. The classical dynamics arises from a 
hyperbolic symplectomorphism on T 2 , which is represented by a matrix A G SL( 2,Z); it can 
be quantized into a family of unitary operators (Uffj^o of finite dimensions ~ hr 1 . One can 
analyze the microlocalization on T 2 of the eigenstates of the operators Uf r For any periodic 
orbit 7 of the classical map, the very special algebraic properties of the operators Un. allows 
to explicitly construct eigenstates with strong scars on 7, with weights w 1 necessarily taking 
values in [0,1/2] [FN04], These special eigenstates were constructed by linear combinations 
of evolved Gaussian wavepaekets localized on 7. 

In the present article we will consider a similar construction for the dynamics provided 
by the geodesic flow on ( M,g ), or more general Hamiltonian flows. By doing so we won’t 
be able to construct eigenstates of P(h), but only approximate eigenstates, which are called 
quasimodes of P(h). Let us now give a precise definition of this notion. 

Definition 1.2. For a given energy level E > 0, we say that a semiclassical family ((pn)n-> 0 
of L 2 normalized states is a family of quasimodes of Pn of central energy E and width f(h), 
if and only if there exists ho > 0 such that 

II (P(h) - E)(p h \\ L 2 {m) < f(h) , Vh G ( 0 , ho] . (1.3) 

A slightly broader definition consists in allowing the central energy to depend on h as well, 
namely considering a sequence {E h ) h with E h —>• E, and requiring 

II (P(h) - E h )y h || L2(M) < f(h) , Vh G (0, ho] . (1.4) 

Remark 1.5. One way to construct quasimodes of center E h and width f{h) is to take linear 
combinations of eigenfunctions of P(h) with eigenvalues in the interval [Eft — f(h), En + f(h)\. 
We will use this trick in the last stage of the proof of Theorem 1.7 below. However, the 
construction of a fully localized quasimode in Prop. 1.9 will proceed differently, namely by 
a time averaging procedure. 

What is the interest of considering quasimodes instead of eigenstates? Except for very 
special systems like the quantized “cat map”, we are unable to give explicit, or even approx¬ 
imate expression of the eigenstates of Anosov systems. On the opposite, we will construct 
explicit quasimodes of certain widths f(h). The larger the width, the less constrained the 
localization properties. Using the fact that the function p(x, f) is the principal symbol of the 
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operator P(h), one can show the following properties of semiclassical measures associated 
with ((p h ): 

• if f{h) = o(l)ft_,.o, then p sc is supported on p~ l (E), and is a probability measure. 

• if f(K) = o(h), then /i, sc must be invariant through the Hamiltonian flow T*. 

For an Anosov flow like the geodesic flow on (M, g), there exist many invariant measures, the 
simplest ones being the Liouvillc measure on p~ l (E), or the delta measures localized along 
a closed orbit 7 C p~ 1 (E). Our Prop. 1.9 will exhibit quasimodes of width Ch/\ log h\ (with 
C > 0 large enough) converging to the measure 5 7 . An easy corollary is that, if the flow 
is Anosov on p~ 1 (E), then for any flow-invariant measure /j on p~ x (E) one may construct a 
sequence of quasimodes of width CK/\ log h\ converging to p. 

In [An06, AN07] it was specified that the semiclassical measures associated with quasi¬ 
modes of width f(h) = ° ( | lci g ft| ) must satisfy the same “half-delocalization” constraints as 
the measures associated with eigenstates. Hence, for such Anosov flows the “logarithmic 
scale” c | lo g seems to be critical for the quasimode width, the constraints on the local¬ 
ization depending on the size of the factor c > 0. In [Brl3] Shimon Brooks studied such 
“logarithmic quasimodes” for the Laplacian on compact hyperbolic surfaces. He showed that 
quasimodes of width c | J) ^ with c > 0 arbitrary small could still be “strongly scarred” on 7: 

Theorem 1.6. [Brl3] Let ( M,g ) be a compact hyperbolic surface, and 7 C S*M a closed 
geodesic. Then for any e > 0, there exists 6(e) >0 and a sequence (<Ph)n-+ 0 of quasimodes of 
center energy E = 1 and width e y lo ^ ^ , converging to a semiclassical measure p sc satisfying 

Tsc( 7 ) > 8(e). 

The methods used in [Brl3] are quite specihc to the setting of compact hyperbolic surfaces. 
Our main objective in this article is to generalize Brooks’s theorem to more general surfaces 
and Hamiltonian flows. Our construction will be mostly local, independent of the global 
dynamics on p~ 1 (E)] essentially, the only assumption we need to make on the dynamics is 
the presence of a hyperbolic closed orbit 7. We will also attempt to optimize the relation 
between the width constant £ > 0 and the weight p sc ({'y})- In principle our methods could 
be extended to higher dimensions, but probably at the expense of this optimization (see the 
Remark 3.5 at the end of section 3). Our main result is the following Theorem. 

Theorem 1.7. Let ( M,g ) be a smooth compact Riemannian surface without boundary. Let 
P(h) € T™(M) be a self-adjoint elliptic pseudodijferential operator, with semiclassical prin¬ 
cipal symbol p(x,f). For a regular energy level E 0 , assume that the Hamiltonian flow 
admits a closed hyperbolic orbit 7 C p~ 1 (E 0 ). 

Choose any e > 0. Then there exists 8(e, 7) > 0, a sequence of energies (E h = E 0 + 
0(h/\ log ft|) and a sequence of quasimodes (Vt)/i^o centered on E h and of width e(h/\ loghj), 
such that any semiclassical measure p sc associated with (fh)h-^o satisfies 

Tsc({ 7» > 8(e, 7). 

Furthermore, we have an explicit estimate for 5 (6,7) in the regime e« 1 : 

*•■”> = ^i7! + 0 ((^) 2 )- (1 ' 8) 

where A 7 > 0 is the expanding rate per unit time along the unstable direction of the orbit 7. 
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In the course of proving this theorem we will first construct logarithmic quasimodes fully 
localized along 7 , that is converging to the semiclassical measure <5 7 , but with a lower bound 
on their width. 

Proposition 1.9. Let (M,g), P(h) and 7 c p _1 (£ , 0 ) be as in the above theorem. 

Then, for any sequence of energies (En = E 0 + O(h))fi ^0 and any constant C 1 > 7 tA 7 , there 
exist a family of quasimodes (h/y e L 2 (M ))^ 0 centered on En and of width C~,(h/\ logh|), 
converging to the semiclassical measure /i. sc = 8 1 . 

Remark 1.10. Notice a difference between the “central energies” En in the Theorem and the 
Proposition: in the latter we are free to choose the energy En in the range Eq + 0(H), while 
in the theorem the energy Ej, is given to us. 

The dynamics covered by our result include various types of dynamics: the Anosov systems 
described above, for which it seems for the moment very difficult to explicitly describe quasi¬ 
modes of widths o(h /\log h \). On the other hand, our result also encompasses 2-dimensional 
Liouvillc-integrablc Hamiltonian flows featuring one hyperbolic orbit: in this case one can 
generally obtain a precise description of quasimodes of width 0(h °°), or even of the eigen¬ 
states of P(h). We will discuss this case a bit in the next subsection. 

1.1. Discussion of related results. The construction of quasimodes which localize along 
closed geodesics, or more generally along invariant submanifolds, has a long history, start¬ 
ing at least with the work of Keller [K58], who initiated the construction of quasimodes 
supported on invariant Lagrangian submanifolds of T*M , using WKB-type Ansatze. The 
methods were further developed during the ’70s and ’80s, at least in the case of completely in- 
tegrable dynamics, by a plethora of authors including Maslov [Mas72], Duistermaat [Duis74], 
Weinstein [Wein74], and Colin de Verdiere [CdV77], through the use of “softer” symplectic 
geometric methods and Fourier integral operators. One could then exhibit quasimodes with 
widths of order 0(h °°), localized on some invariant Liouville-Arnold torus. 

Babich and Lazutkin [BabLaz 68 ] constructed 0(h°°) quasimodes concentrating along a 
single closed geodesic of elliptic type (that is, the linearized Poincare map of the orbit 7 is an 
elliptic matrix). The works of Voros [Vor76], Guillcmin-Weinstein [GuiWein76], and Ralston 
[Ral76] each carried out this idea via different methods; in particular Ralston constructed 
maximally localized eigenstates in the form of Gaussian beams , which are minimal uncer¬ 
tainty Gaussian wavepackets transversally to the orbit. These constructions were specific to 
the case of elliptic orbits, albeit with nonresonance conditions. 

In the case of a hyperbolic orbit, the construction of quasimodes using Gaussian 

wavepackets usually breaks down, due to the fast spreading of the wavepacket along the 
unstable direction. However, for certain Liouville-integrable systems, the particular structure 
of the unstable and stable manifolds of 7 allows to construct 0(h°°) quasimodes as WKB 
states along the various branches of these manifolds. This construction, peformed by Colin 
de Verdiere and Parisse [CdVP94] in the case of a surface of revolution, shows that some 
family of eigenstates of the Laplacian converge to the semiclassical measure <5 7 , but the 
convergence is rather slow: the weight of the measure | <pft(x)\ 2 dx on the complement of a 
small neighbourhood of 7 decays at the rate | log hi” 1 . We will come back to the Colin de 
Verdiere-Parisse quasimode construction in subsection 5.3. Toth studied specific integrablc 
systems in higher dimension [To96, To99], and showed that some families of eigenstates may 
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converge to 5 7 for 7 a hyperbolic orbit. The inverse logarithmic decay of the measure away 
from 7 was shown to hold in greater generality by Burq-Zworski [BZ04]: their Theorem 2’ 
implies that in any dimension, a quasimode of width ch/\ log H\ with c > 0 small enough must 
have a weight > C -1 | log h\~ l outside a small neighbourhood of a hyperbolic closed orbit 
7. These results were generalized to semi-hyperbolic orbits by Christianson [Chr07, Chi'll], 
including the case of a manifold with boundary. 

In the Quantum Chaos physics literature several studies were devoted to constructing 
various forms of localized states, which could hopefully mimick the “scars” numerically ob¬ 
served on certain eigenfunctions, or at least exhibit significant overlaps with these eigenstates. 
These localized states are generally based on Gaussian wavepackets localized on a point of 
7. De Polavieja et al. [dPBB94], and then Kaplan and Heller [KH99] used time averaging 
of such a wavepacket to construct a localized quasimode, but neither sets of authors tried to 
optimize the time of evolution or the width. On the opposite, a series of works by Vergini and 
collaborators constructed so-called “scar functions” localized along 7. They first considered 
a Gaussian beam along 7 [VCOO], then used several procedures to improve the width of the 
quasimode, keeping the latter localized near 7 [VC01, VS05]. In particular, in [VS05] they 
improved the Gaussian beam through a time averaging procedure, and optimized the width 
by evolving up to the Ehrenfest time and by carefully selecting the weight function. 

In the mathematics literature, Brooks [Brl3] also constructs localized logarithmic quasi¬ 
modes in the course of proving his Thm 1.6; his method also relies on time averaging, but 
starts from “radial states” specific to the geometry of hyperbolic surfaces, instead of Gauss¬ 
ian wavepackets. As Brooks did not try to optimize his constants, one of our aims in this 
article is to make those precise. 

1.2. Outline of the proof. Most of the article deals with the proof of Prop. 1.9, namely 
the construction of the fully localized quasimodes. Our proof is essentially identical to the 
procedure used by Vergini-Schneider in [VS05]. The main difference is that the semiclassical 
tools we use (e.g. the Quantum Normal Form) are mathematically rigorous, including our 
remainder estimates. The extension to the partially localized quasimodes of Thm 1.7 is 
similar in spirit to that of [Brl3], albeit with a better control of the constants. 

We now present an outline of the proof of Theorem 1.7, which starts by constructing the 
fully localized logarithmic quasimodes of Prop. 1.9. 

1) Quantum Normal Form The first task is to choose convenient canonical coordinates 
near the orbit 7, namely normal form coordinates, which map the neighbourhood of 7 = 
7 (E 0 ) C p - 1 (E 0 ) into a neighbourhood of the circle {r = x = £ = 0, s G T} inside the model 
phase space T*(T S x M x ) spanned by the coordinates (s, x ; t , £). In these coordinates, the 
Hamiltonian p takes the following form, at the order 2 N in the transverse coordinates: 

p(s,x]T,f) = f(r) + \(t)x£+ Y Qa(r){x^) a + 0((x ,£) 2N+1 ), 

2 <ct<N 

where /(r) = E 0 + -fr + (9(r 2 ), with T 0 the period of 7 = ^{Ef) and A(0) > 0 the expansion 
rate along the unstable direction. This normal form can be extended to the quantum level, 
following works of Gerard-Sjostrand and Sjostrand [GS87, Sj 02 ]. The result is that, microlo- 
cally near 7, the operator P(h) is unitarily equivalent to a Quantum Normal Form (QNF) 
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operator P(h), a pseudodifferential operator on I x I of “simple” form; in particular, its 
principal symbol of P(h) is given by p above. 

In this outline we will rather consider a simplified normal form to focus on the important 
ideas. Namely, we will take the following quadratic Hamiltonian on the phase space T*(T x 


P 2 (s,x;t,£) = Eo + ^ + Xx^,, 

with A = A(0). The stable (resp. unstable) subspace corresponds to the £-a.xis (resp. x-axis). 
The corresponding quantum operator is obtained through the Weyl quantization of p 2 on 
T x I 9 (s, x): 

P 2 (h) = Op^(p 2 ) = E 0 + —hD s + A (xhD x — ih/ 2), s G T = R/Z, igK. 

To 

On the right hand side we recognize the selfadjoint ID dilation operator (xhD x — ih/ 2). 

2) Propagation of transverse Gaussians: The state with which we start has the shape 
of a Gaussian beam: it is a minimal uncertainty Gaussian wavepacket transversely to 7 , and 
a plane wave along the direction of 7. In our normal form coordinates (s,x), this state has 
the form 


ipo (s,x) 


(^yi /4 ex p(2fvrns) exp 



for some choice of n G Z. 


(1.11) 


This state has, with respect to the model Hamiltonian P 2 (h), an average energy 0 , P 2 {h)ipo) = 
(ipo, -F^o) = E 0 + ^r-h, and an energy width 0(K). We hrst propagate this Gaussian beam 

using the Schrodinger equation; for the model operator P 2 (h) the propagation is explicit, 
and preserves the factorization of variables: 

AM = [e- w/i «(s,!) = e-"*' 1 pp ( - jsj) ■ (1.12) 

The longitudinal dependence e 2mns of the state is unchanged, but its transverse part is now 
a Gaussian of spatial width ~ e tA h 1,/2 , which grows exponentially when t > 0 (recall that the 
unstable manifold is the x-axis). The state remains localized in a microscopic neighbourhood 
of 7 until the Ehrenfest time T\ = d-d log H\. Notice that, with respect to the full geodesic 
flow on S*M , T\ is only the local Ehrenfest time associated with the orbit 7. 

To control the evolution of ipo with the full QNF Hamiltonian P(h), we will use the results 
of Combescure and Robert [CR97] who give an explicit procedure to calculate e -* tp (k> / h ip n in 
terms of an asymptotic expansion in powers of h (see also [Hag80]). Although their expansion 
usually breaks down at the time ~ \T\, the special structure of the QNF P(h) allows here 
to keep the evolution under control until T\. Besides, up to this time the expansion is 
dominated by the state i/t of (1.12). 

3) Time averaging: Obviously, each evolved state i/> t is a quasimode of P 2 with the 
same center and width 0(h) as ipQ. On the other hand, one can decrease the width of 
the quasimode by linearly combining these evolved states, namely averaging over the time. 
Following the approach in [VS05], we consider the Ansatz: 

/ OO 

Xit/T) 

-OO 


(1.13) 
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where a £ C^°((—1,1)) is a smooth cutoff, and T > 0 is a (large) time. Both a and T will 
be optimized later. To compute the width of this quasimode, we must compare the norms 
of the states (P — E h ) x V h and Both computations rely on first estimating the overlaps 
(i/y, tV) between evolved states. For our model Hamiltonian P 2 the overlap is explicit: 




, 1 For the evolution through the full QNF operator P(h), we will see 

W 2 cosh((t— t')X) 


that these overlaps decay essentially in the same way. An integration by parts shows that 


(P - E h ) [ a (t/T) e itE ^ h Adt='^ [ X '(t/T) e itE ^ h ^ dt . (1.14) 

Jr 1 Jr 


Hence, the larger the time T, the smaller the prefactor in the RHS. Because the evolved 
states remain localized until the time T\, we will choose the cutoff time T = (1 — e')T\ 
for some small s' > 0. The prefactor of the above right hand side explains the order of 
magnitude C of the width. To minimize the constant C we also need to optimize the 
cutoff function a- This optimization completes the construction of fully localized quasimodes 
in Prop. 1.9. 

4) Projections onto logarithmic intervals: After constructing our logarithmic quasi¬ 
mode Tft with center Em width C' 7 , and which is microlocalized on 7 , we project T j, on 
an energy interval of width 2eh/\ log h\ centered near Em this projection is then automati¬ 
cally a quasimode of width eh/ | log Ti\. Elementary linear algebraic arguments show that, if 
the center E' n of this interval is carefully chosen, the projected state is nontrivial, and has a 
positive weight on a small neighbourhood of 7. 


Remark 1.15. The compactness of M has not been used in the proof of Prop. 1.9, but only 
in the final step of the proof of the Theorem, where we need the spectrum of P{h) to be 
discrete. This discrete spectrum is also the reason why we required P(h) to be semiclassically 
elliptic. Actually we need this discrete spectrum only in a neighbourhood of the energy 
E 0 , so the result of Thm 1.7 can be generalized, for instance, to the case of Hamiltonians 
p(x,£) G S m (T*M 2 ) such that the energy layer p _ 1 ([P 0 — e, E 0 + e]) is compact. 

Acknowldegements. The work in this article began while S.E. was a resident at the Insti- 
tut de Hautes Etudes Scientihques, and was completed while being supported by an AMS- 
Simons Travel Grant. S.N. was partially supported by the grant Gerasic-ANR-13-BS01- 
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2 . Semiclassical preliminaries 

In this sections we recall the concepts and definitions from semiclassical Analysis we will 
need. The notations are drawn from the monographies [DimSj99, Z12], 

2.1. Pseudo differential operators on a manifold. Recall that we define on M? d the 

following class of symbols for m G R 2 : 

S m (R 2d ) = {a G C°°(R 2d x (0,1]) : |<9"dfa| < C a ^) m ~ m }. (2.1) 

Symbols in this class can be quantized through the h-Weyl quantization into the following 
pseudodifferential operators acting on u G iS(M d ): 

Op«(°)'“W= / 2 /n) d / 


( 2 . 2 ) 
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One can adapt this quantization procedure to the case of the phase space T*M, where M 
is a smooth compact manifold of dimension d (without boundary). Consider a smooth atlas 
(fhVi)i=i,...,L °f M, where each /) is a smooth diffeomorphism from VJ C M to a bounded 
open set Wi C R d . To each fi correspond a pullback /* : C°°(Wi) —> C°°(Vi) and a symplectic 
diffeomorphism fi from T*Vi to T*IT): 

fi ■ (z,0 ^ (Dfi{x)~ 1 ) T t > ) . 

Consider now a smooth partition of unity (<pi) adapted to the previous atlas (fi,Vi). That 
means Yh <Pi — 1 and <pi G C°° (V)). Then, any observable a in C°°(T*M) can be decomposed 
as: a = a/, where a; = a<^/. Each cq belongs to C°°(T*Vi ) and can be pushed to a function 

di = ( f l ~ 1 )*ai G C' 00 (T*W / i). We may now define the class of symbols of order m on T*M 
(after slightly abusing notation and treating (x,£) as coordinates on T*Wi) 

S m (T*M) = {a G C°°(T*M x (0,1]) : a = a h such that (2.3) 

i 

~ ai G S m (R 2d ) for each /}. 

This class is independent of the choice of atlas or smooth partition. For any a G S m (T*M), 
one can associate to each component di G S m (M? d ) its Weyl quantization Op^(a/), which 
acts on functions on M 2cZ . To get back to operators acting on M, we consider smooth cutoffs 
ifji G C^°(Vi) such that tpi — 1 close to the support of ipi, and dehne the operator: 

Op»(o)« = x (/.-OpJ-fSiK/r 1 ) - ) (* x u), ii e C“(M). (2.4) 

i 

This quantization procedure maps (modulo smoothing operators with seminorms 0(h°°)) 
symbols a G S m (T*M ) onto the space \1/™(M) of semiclassical pseudodifferential opera¬ 
tors of order m. The dependence in the cutoffs tpi and ^ only appears at order 
([Ho85, Thm 18.1.17] or [Z12, Thm 9.10]), so that the principal symbol map <t 0 : 

S m (T*M)/hS m ~ 1 (T*M) is intrinsically defined. Most of the rules and microlocal properties 
(for example the composition of operators, the Egorov and Calderon-Vaillancourt Theorems) 
that hold on W l can be extended to the manifold case. 

An important example of a pseudodifferential operator is the semiclassical Laplace-Beltrami 
operator P(h) = In local coordinates (x;£) on T*M, the operator can be expressed 

as °Ph(\C\l+H'E j b j( x )£j+ c (x))+h 2 d(x)) for some functions bj , c, d on M. In particular, its 
semiclassical principal symbol is the function |£| 2 G S 2 (T*M). Similarly, the principal sym¬ 
bol of the Schrodinger operator -yA s + h(x) (with V G C°°(M)) is |£| 2 + V(x) G S 2 {T*M). 

We will need to consider a slightly more general class of symbols than the class (2.1). 
Following [DimSj99], for any 0<<5<l/2we introduce the symbol class 

S™(R 2d ) d = {a G C°°(R 2d x (0,1]) : \d^a\ < 

These symbols are allowed to oscillate more strongly when h —> 0. All the previous remarks 
regarding the case of 5 = 0 transfer over in a straightforward manner. This slightly “exotic” 
class of symbols can be adapted on T*M as well. For more details, see [DG13, Section 3] or 
[Z12, Section 14.2], 
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3. Normal form around a hyperbolic trajectory 

Although the 2-dimensional quantum Birkhoff normal form (abbreviated by QNF) that we 
present below from [Sj02] was proved in the setting of real analytic Hamiltonians on T*R 2 , 
its proof directly transfers into the C°° setting, and we will present a sketch of it for the sake 
of the reader. 

Let us recall our setting. (M, g ) is a smooth surface without boundary and P(h) is 
a pseudodifferential operator on M with principal symbol p € S m (T*M) assumed to be 
independent of h, and formally selfadjoint on L 2 (M, dg ). We assume that for a regular energy 
level E 0 (meaning that dp does not vanish on the energy shell p~ 1 (E 0 )), the Hamiltonian 
flow exp{tH p ) admits a hyperbolic periodic orbit 7 = 'y(Eo) C p _1 (-E 0 ), with period T(E 0 ). 
The hyperbolicity implies that this orbit is isolated in p~ l (E 0), and that it belongs to a 
smooth family {7(TJ)} of hyperbolic orbits, with which each inside of p~ 1 (E). The linearized 
Poincare map of 7 has eigenvalues A (E 0 ), A(Ao) _1 with |A(_E 0 )| > 1; A (E 0 ) is positive (resp. 
negative) if the local unstable manifold of 7 is orientable (resp. nonorientable). In both cases 
we denote |A(i? 0 )| = exp(T(E 0 )\( y E 0 )), so that A (E 0 ) is the expansion rate per unit time 
along the unstable direction. This parameter will play an important role in the following. 

Denoting by a the Liouville 1-form on T*M, which reads a = Yl?j=i V.j dlJj i n local coordi¬ 
nates ( ip, Pi ) it. We call 

p = <p(E 0 ) = / a, the action of the orbit 7 (E 0 ). 

J 7 ( E 0 ) 

Proposition 3.1. [Sj02] i) There exists a Birkhoff normal form around the orbit 7(^0). 
Namely, for any integer N > 0, there exist local symplectic coordinates (s,r,x,f) G T*(S i x 
R) and a smooth local canonical transformation 

n = n N : neigh(70 = {(s, 0, 0, 0), s G S' 1 }, T *(T x R)) —>• neigh(7(E 0 ), T*M) , 

single-valued if the unstable manifold of'y(Eo) is orientable, and otherwise double-valued with 
k,(s — 1, r, x, f) = k(s, t, —x, —£), such that 

p N = pon = f(T)+\(T)xf+q iN) (r,xf)+0{(x,f) N+1 ), f N \ r ,xf)= ^ q a (r)(xf) a . 

4<2cKiV 

(3.2) 

Here f(r) = E 0 + fr + 0(t 2 ), and to each value t e [—e, e] corresponds a periodic orbit 
of energy Efr), with E( 0) = E 0 . By a slight abuse of notation we denote A(r) = A (E(t)) 
the unstable expansion rate for the orbit 7 (E). The remainder means that \0((x, A 

C(\x\ + |^|) Ar+1 in the indicated neighborhood ofT *(T x R). 

ii) Recall that tp is the action of the orbit 'y(Eo). Given h e (0,1], let S°^™(T x R) 
be the space of smooth functions u(t,x ) on R x R, in the Schwartz class in the second 
variable, with the periodicity property u(s + l,x) = e~ lip ^ h u{s,x) in the orientable case, resp. 
u(s + l,x) — e~ lip ^ h u{s, —x) in the non-orientable case. 

There exists a corresponding quantum. Birkhoff normal form, namely a semiclassical Fourier 
integral operator Un : «S“^°(T x R) —* C°°(M) quantizing microlocally near 70, microlo- 
cally unitary near'y(Eo) X70 (using the natural L 2 (TxR) structure). This operator conjugates 
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P = P(h) to the following Quantum Normal Form operator: 

U* N PU N = P^ + R N+1 . (3.3) 

Here P {N \ R N+1 are pseudodifferential operators on S^/niY x M), microlocally supported 
near 7 0; with Weyl symbols 1 

N 

P (N) (r,x,£;h) = £/^(r,*,0, r N+1 (s,T,x,£-,h) = 0((h,x,f) N+1 ) , (3.4) 

3=0 

where p^ = p n above, and the higher order terms are of the same form. Moreover, 
p {N+1) ( T , x, Q h) - p^ N) ( t, x, Q h) = 0((h,x,f) N+1 ). 

Proof. (Sketch, after [GS87, Sj02]) 

i) We will only consider the orientablc case. The first step of the proof lies in the con¬ 
struction of suitable symplectic coordinates kn(s, t; x, £) = (y,r)) such that the classical 
Hamiltonian p = <Jo (P) in these coordinates reads as in (3.2). That is, we construct the 
classical normal form near the orbit 'y(Eo). Let us recall this construction. 

As mentioned above, the hyperbolicity of 'j(Eq) implies that in some range E G [Eq — 

e 0 ,E 0 + e 0 ], the flow admits a hyperbolic orbit 'y(E), so that T = f U|B-_Eo|<eo forms 
a smooth symplectic submanifold. If we call F±(E) the local unstable/stable manifolds of 
7 (E), which are smooth immersed Lagrangian leaves, their unions T± = U|E-E 0 |<e 0 T±(-E) 
are smooth involutive outgoing/incoming submanifolds of codimension 1. One can use 
the implicit function theorem to define a smooth function £ near 7 (Eq) such that locally 
T+ = {£ = 0}, and a symplectically conjugate coordinate x with x [r_= 0. We may 
also setup canonical coordinates (s, r) on T, with s multivalued, such that each orbit 7 (E) 
corresponds to a circle {r = r(E), s G [0,1)}. This set of coordinates on T can be ex¬ 
tended to a neighbourhood of T, such as to produce a system of canonical coordinates 
(y,rj) = h 2 (s,t, x, £) with r = r(E) constant on each submanifold T ±(E). In these coordi¬ 
nates near 7(i7o), the Hamiltonian is of the following form: 

P 0 «2 (t, t; x, £) = /(r) + A(t)t£ + 0((x, £) 3 ) , 

where the function r hg f(j) is the inverse function of A ha t(E). In order to normalize the 
higher order terms in (x,£), we can iteratively construct functions (Gj(t,r,x,£))j> 3 (which 
are homogeneous polynomials of degree j in (x, £), satisfying certain solvable ODEs), and 
use the Hamiltonian flows e tH °i they generate, such as to obtain, at the order N, the normal 
form: 

p N = po K 2 o e H °3 O • • • O e Ha N = /(r) + A(r)x£ + ^ g a (r)(x£) JV + 0((x, £) JV+1 ), 

2 <a<N 

which is the desired Birkhoff normal form. The induced change of coordinates kn = K2°e HG 3 o 
... o e H °N maps a neighbourhood of 70 G T*( Tx M) to a neighbourhood of 7(^0) G T*M. 
Notice that kn and kn+i have the same Taylor expansion in the coordinates (x, £) up to 
order N. 

^hince T = R/Z, the Weyl quantization is intrinsically defined on T*(T x R) by pullback from the 
quantization on T* R 2 . 
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ii) One can adapt this normal form construction to the quantum framework (a general dis¬ 
cussion on quantum normal forms can be found in [Z12, Chap.12]). The task is to construct 
a semiclassical Fourier Integral Operator (FIO) Uq quantizing the canonical transformation 
Kn. The construction of the phase function of U$ shows that this operator maps the space 
C' 0O (M) to the space x M) of functions which are Z-periodic functions in the variable 

s, twisted by the phase e llf ^ h (see for instance [GS87, Eq.(2.9)]). The FIO U 0 is elliptic from 
a neighbourhood of q 0 to a neighbourhood of 7 (E 0 ), and can be constructed such as to be 
microlocally unitary between two such neighbourhoods. 

Conjugating our quantum Hamiltonian P(K) with Uq, we obtain a pseudodifferential oper¬ 
ator P 0 (h) l = UqP{K)Uq on T x M), whose principal symbol is of the form Pn(s, t,x,£). 
The subprincipal symbol of Pq ( K) is, a priori, an arbitrary function of (s, t;x, £). By itera¬ 
tively solving a sequence of transport equations, we may correct the symbol of the FIO Uq 
into a FIO Un, so that the symbol of the conjugated operator Pn(K) = U^P(H)U]y takes the 
form (3.4). 

Our main point is that the above construction, originally presented in the case of analytic 
symbols and operators on the Euclidean space [GS87, Sj02], can be generalized to the case 
of smooth objects on a smooth manifold (the FIO Uq can be constructed with a real valued 
phase function). Since the construction is local near 7(E 0 ), topological properties of the 
manifold M (e.g. the homotopy class of 7 (E 0 )) do not come into play. □ 

Remark 3.5. Similar quantum normal forms were established in higher dimension. The refer¬ 
ences [Gui96, Zcl98] are specific to the homogeneous microlocal setting. A higher dimensional 
formulation of the QNF we use was given in [GuiPaul09], albeit the proof for the case of a 
hyperbolic trajectory is only sketched. A similar QNF was also derived in [Chr07], but at 
the cost of a remaining elliptic factor depending on the longitudinal variable. 

Using these higher dimensional QNF, one can probably generalize the quasimode con¬ 
struction we are presenting below. The main extra difficulty is to compute the corresponding 
square norm of the quasimode T as performed in Lemma 5.4. This computation requires 
to understand well the overlaps between evolved coherent states in higher dimension (the 
generalization of (5.7)), which will depend on the spectrum and Jordan structure of the 
linearized Poincare map of 'y(Eo). Since obtaining sharp constants is one objective of the 
present article, we restrict ourselves to the 2-dimensional case. 


4. Propagation of a Gaussian wavepacket at the hyperbolic fixed point 


The construction of our quasimode will be performed on the Quantum Normal Form side, 
that is on the model space T x K. As explained in §1.2, it will be based on an initial 
state which is a plane wave in the longitudinal direction, and Gaussian along the transverse 
direction. We must nevertheless take into account that the QNF operator acts on twisted- 
periodic functions. From the action (p(E 0 ) and h G (0, 1] we setup the number 


def 

Vh = 



<p{Eo) 

27i h 


1 


where [s] represents the integral part of s G R. . 


It satisfies e liph = e l ‘ p ^ E °' > l h , so the spaces of twisted periodic functions can be denoted by 
S°J h n °{S l x M). With this twist being taken into account, our initial state will be of the same 
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form as in (1.11): 

ipo{s,x) = e 


4 ?! i(2nn—(pti)s 


7T 


ny / 4 


exp 


x~ 

(-for some arbitrary n £ Z. (4.1) 

Z/t/ 


Due to the parity of a; ha- e~ x2 ^ 2h , this function belongs to both spaces S^ l °(T x R), and it 
will allow us to treat both the orientable and nonorientable cases. 

We will select the index n to be uniformly bounded when ft —> 0, so that 


OPft(/( r ))^o = f(hD s )y >o = f(ft(2irn - pn))y o satishes f(ft(2nn - (p h )) = E 0 + 0(ft). 

We will evolve the state y 0 through the Schrodinger equation generated by the QNF operator 
P (A 4 described in Prop. our task will now be to describe the states 


= e * p(JV) *A^ 0 , for times \t\ < C \ logh| . 
The Weyl symbol p ^ of can be written 


N 

P iN \r, x, £-h) = Y^ Q a (P K) , 

a=0 


(4.2) 


where the symbols q a (r, ft) expand as 

N 

q a {»] ft) = ^ ft 1 g“(«) ? with the special values q°(r) = /(t), ql(r) = A(r). 

i =0 

The major advantage of this QNF is the separation between the variables (s, r) and (x, £). 

As a result, the Schrodinger evolution can be reduced to a family of 1-dimensional problems. 

Indeed, since e l ^ 2wn ^ v>h ^ s j s an eigenfunction of ftD s , the state ip t can be factorized into 
longitudinal and transversal parts, 

y>t N \s,x) = e~ ltq °^ h e l ^ 2wn ~ (ph ^ s (p[ N \x) , where we use the shorthand notation q° = q°((2nn—ipfi)ft', ft) ■ 
Notice that q° = E 0 + O(ft). The transversal part ip[ N ^ satishes the ID Schrodinger equation 

m,4 m = Q ,m (K) Q im = Op 

N 

q {N \x, f; ft) = ^ Q a j. wliere we took (f = <?“((2vm - (p h )ft ; ft) . (4.3) 

a=l 


We thus end up with analyzing the evolution of the ID Gaussian state tpo (x) = (t/I) -1 / 4 e - ^ 
under the effective ID Schrodinger equation. 

Let us truncate the symbol q (A 4 to the quadratic order, that is keep from (4.3) the term 

g(A) — q 1 x £ . where q 1 = q l {( 2 -nn — (pa)ft ; ft) = A(0) + Oift). (4.4) 

The evolved state through this quadratic operator is easily expressed: 

exp ( - l -^Q { q N) )(fo = PtqiPo , (4.5) 

and we used the unitary dilation operator T>g : L 2 (R) -A L 2 (R): 

Vpu(x) = f exp (—i/3 Op %(x£)/ft)u(x) = e~^^ 2 u{eT^x). 


(4.6) 
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The state V tq up 0 is known in the literature as a squeezed coherent state [CR97]: it is a 
Gaussian state, but with a different width from that of ipQ. This state is still centered at 
the origin in phase space, (x, £) = (0,0), which is a fixed point for the classical evolution 
generated by the Hamiltonian q^ N P 

For N > 2, the operator includes a nonquadratic part 

N 

C = s™ - 9 (1) = £ (4-7) 

a=2 

We show below that, in the semiclassical limit, this nonquadratic component of Q ^ will 
induce small corrections to the state (4.5) evolved through the quadratic part. To justify 
this fact we will use a Dyson expansion, like in the general treatment of [CR97, Section 
3], which considered the evolution of coherent states through arbitrary Hamiltonians. Here, 
the special form of q( N l will facilitate our task and will produce smaller remainder terms. 
Remember that the state (4.5) is a squeezed coherent state. The corrections due to the 
nonquadratic part are linear combinations of squeezed excited states. 


Squeezed excited states The initial coherent state (po is the ground state of the stan¬ 
dard ID Harmonic oscillator ( hD x ) 2 + x 2 . We will call (<p m ) m >i the m-th excited states, 

which are obtained by iteratively applying to (po the “raising operator” a* = f OPft(^=|) and 
normalizing: 


*■Pm 


<Po 


<Pr 


(X = 


(irh) 1 / 4 2 m / 2 \fm\ 


H m (x/h 1/2 ) e 


2h 


(4.8) 


where H m (•) is the m-th Hermite polynomial. By applying the unitary dilation operator Vp, 
we obtain a family of squeezed excited states {T>p(p m ) m >i. 

Following the strategy of [CR97], we will show the following approximate expansion for 
e ~itQ (N '>/h)tp 0 j n terms of squeezed excited states. 


Proposition 4.9. For every l e N, there exists a constant Ci > 0 and time dependent 
coefficients c m (t,h ) G C for 0 < m < 21, such that the following estimate holds for any 

he ( 0 , 1 ]; 

21 

Vt e R, ||e^ W /Vo - Vtqupo - Y, Fn(t, h)v tq up 2m \\ < a (\t\h ) l+1 . (4.10) 

m =0 

The coefficients c m (t,h ) are polynomials in ( t,h ), with degree at most l in the variable t. 
Besides, co(t, h) = O t (h), c\(t, h) = O t {h 2 ), c m (t , h) = O t (h^ m+1 ^ 2 ^) for 2 < m < 21. 

Remark 4.11. The result presented in [CR97, Thm 3.1] concerned the evolution of tp o w.r.t. 
an arbitrary quantum Hamiltonian; in this general case the remainder usually grows ex¬ 
ponentially with the time. The present polynomial growth of the remainder relies on the 
normal form structure of Q^ N \ namely the fact that (fisp is a sum of powers of Op^ (^)i 
which commutes with the quadratic evolution. 

Proof. Like in [CR97], we want to compare the full evolution = U(t)ipo = e~ ltQ[N) ^ h (po 
with the quadratic one, U q (t)ipo = e~ lt ® q ^ h <po- The comparison is based on the Dyson 
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expansion 

u(t) - U q (t) = ±-J*U(t- U) QW u q (h) dt \. (4.12) 

This is the first order Dyson expansion (case 1 = 0). Let us show that QlPP commutes with 
Q[ N P Indeed, each operator Op^((:r£) a ) is a linear combination of powers of Op^ 

Lemma 4.13. For any a > 2, there exists absolute constants {c a ^)i<k<[a/ 2 } such that 

[a/2] 

<«(«)“) = (OpJfLOr + E c«.A'" 2 ‘(OpJ'(x{))"- 2 ‘ . (4.14) 

fc=l 

This Lemma can be proved by induction, using the fact that the ^-expansion of the Moyal 
product (x£) terminates at the order h 2 . As a result, commutes with and 

therefore with U q (t x )\ 

U(t ) - U q (t) = ±.Ju(t~ h) u q (h) gj dh . (4.15) 

We want to apply this operator to (po. We can represent Qn^Vo as a linear combination of 
excited coherent states. Indeed, using the expression Op n {xC,) = ih((a*) 2 — a 2 ) in terms of 
the raising and lowering operators a*, a, we may write: 

[a/2] 

OpHK) 0 ) = h a (i a ((a *) 2 - a 2 r + YP i a ~ 2 k Ca,k((ci *) 2 - a 2 ) a ~ 2k 

k= 1 

Using the commutation relation [a, a*] = 1, the RHS can be rewritten in “normal ordering”, 
that is as a linear combination of terms (a*) 2/3 a 27 , with 0 < /3,7 < a. Since tpo satisfies 
acpo = 0, we only keep the terms with 7 = 0. Using the definition (4.8) of the excited 
coherent states (p m , we get the following expression: 

Lemma 4.16. For any a >2, there exists coefficients {d a ,k, 0 < k < [a/2]} such that 

[a/2] 

Op R ((a;f )°Vo = fr a Y da ’ k F2a-4k ■ (4.17) 

k= 0 

From (4.7), the state Qn^Vo is a linear combination of 0 < m < N}. By inspection, 

we see that the coefficients in front of (p 2m have the following orders in h: 

QnqF o = ^(^)( < 7’0 + ^ 4 ) + 0(h 3 )(ip 2 + <Pe) + + • • • + 0(h N )ip 2 N ■ 

From (4.5), the action of U q (t \) results in a multiplication by e~ ltiq °/ h , and the replacement 
of p >2 m by the squeezed states V tiq up 2 m ■ Taking into account the factor i/h in front of the 
integral (4.15) and the unitary of U(t), we get the simple estimate 

\\(U(t)-U q (t))<p 0 \\ Lm <Cth, Vhe ( 0 , 1 ], V/eM. 
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In order to improve the description of U(t)(po, we shall expand U(t) — U q (t ) into a Dyson 
expansion of higher order. For any order l > 1, this expansion reads: 

l -t 

= e pmt / ■■■ ■ --Qiyu^t o <#, .. .*, 

j=1 \ in ) Jo Jt i Jtj-i 

+ /'/'••• r d * - < 1 + 1 ) 0 !,"’^,(* 1+1 - miy ■ ■ ■ Q ( ,yu 2 (to *,... *, +1 

Using the commutativity of with U q (s), this simplifies to 

c/(«) - £/,(«> = e u MQiyy + tUst J‘ j,mt ~ t l+1 ) u q (t,+o (Qxyy +i *,+,. 

(4.18) 

When applied to <^ 0 , each j-terrn on the RHS leads to a linear combination of squeezed 
excited states 0 < m < Nj}, with the following estimates on the coefficients: 

fi 0 = —4 + ' ' ‘ + <^ 0 ) + 0(^ +1 ) (<^4j+2 + ^Aj-2 + ' ' ' + <^ 2 ) + 

+ 0(h J+2 )994 :; ' + 4 + • • • + 0(h^ N 

Hence, each j-terrn in the equation (4.18) has a norm 


where the implicit constant only depends on the coefficients q a of Q^ N \ For the same reasons, 
the remainder term in (4.18), when applied to (po, has a norm bounded by 

\-J^f |t/« - t I+ i) U,(t w ) (O'+Vo*,+‘1 < C, (|(|ft)«. (4.19) 

In the ji-surn, each squeezed excited state V tq iLp 2m has a coefficient given by a certain poly¬ 
nomial We won’t need to analyze this polynomial in detail, but only give partial 

information. Because we will deal only with times |f| < C\ logh|, the size of each polynomial 
c m (t, h) will be guided (up to a logarithmic factor) by the term with the smallest h-power. 
Beyond the principal term U q {t)(po, the Dyson expansion will be of the form 

'Dtq 1 \O t (h)(<p 0 + 974) + O t (h 2 )((p 2 + </?6 + ^ 8 ) + ^t(h 3 )(</ 7 10 + ^12) + O t (h 4 )((pi4 + <Pi§) + . . . ] , 

where the coefficient in front of each (p m is a polynomial in t which vanishes when t = 0. In 
the range 3 < j <1, the terms are of the form O t {pj)((p 2 ( 2 j-i) +</? 4 j). This ends the proof of 
the Proposition. □ 


After tackling the ID evolution, we can now reconstruct the full evolved state 
4 N) (s,x) d = e~ itp{N)/h Ms,x) = e- itq ° /h e i{2 ™-^ s <p { t N) (x ). 


(4.20) 
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4.1. Microlocal support of the evolved state. The expansion in Prop. 4.9 shows that 
the state is under control until polynomial times |i| ~ fC 1 . However, for such large 
times the coherent states T>tq 1 T 2 m will be very delocalized. Because the normal form is valid 
only in a small neighbourhood of 70 and it involves remainders 0((x, £) N+1 ), it is important 
to keep our states to be microlocalized in a microscopic neighbourhood of y 0 . For this reason, 
we are forced to bound the time of evolution in a precise logarithmic window. Namely, we 
select a small s' G (0,1), and for any h G (0,1/2] we define the local Ehrenfest time 

« (l-eQUogfil (4.21) 

2A v ; 

where as above A = A(0) is the expansion rate of 70 . 

Proposition 4.22. Take 0 G C'/°(T*(T x M)) with 0 = 1 in a neighbourhood of 70 , and 

denote its rescaling by © a (s, r; x, f) = f 0(s, r/a; x/a, £/a). Then, for any power M > 0, 
there exists Cm > 0 such that 

ll[ 0 p»(e». 7 .)-/]d W, ||y <C M h M , fie ( 0 , 1 ], (4.23) 

uniformly for times t G [— T e /, T e >]. 

The above property could be abbreviated as ||[Op R (0 Re '/ 3 ) — h]fy N} \\m = 0(h °°), uni¬ 
formly in the time interval. Shortly speaking, the state fy N} is microlocalized in any If' 
neighbourhood of 70 

Proof. We will check that all the squeezed excited states T > t q 1 T 2 m, 0 < m < M , are microlo¬ 
calized in the same neighbourhood of (0, 0) G T*K, for t in this time interval. The state 
'Dtqi ( p m (x) is a Gaussian of width e q t h 1 ^ 2 , decorated by a polynomial factor. For |i| < T e t, 
this width takes values: 

e qlt h 1/2 < e qlT *'h 1/2 = e C+OW)T £ , h l /2 = h e'/2+o(H ) ; 

hence it remains microscopic. Consider a cutoff x £ [0,1]) supported in [—2, 2], equal 

to unity in [—1,1], and define Xa( x ) l = xi x / a )- For a > h £,//3 , a direct estimate of the 
Gaussian integral shows that for any M > 0, 

\\( Xa -i)v tq ^ m \\ L2 = o(h M ), 

uniformly in the time window. This shows that T> tq up m is microlocalized inside the strip 
{|a:| < h £ '/ 3 } C T*M. The semiclassical Fourier transform leaves the states ip m invariant 
(up to a constant factor), and inverts the dilation operator: VpifFf'a) = TrffD^yv). As a 
result, the above computation shows that T> tq i(p m is also microlocalized inside the horizontal 
strip {|£| < h £ '/ 3 }, uniformly for \t\ < T E >. These position and momentum mi crolocalizations 
imply that V tq up m is microlocalized inside the square {|x|, |^| < h £ Z 3 }. Hence, for any 

6 G C/°([—2,2] 2 ) with 9 = 1 in [—1,1], rescaled into 9 a (x,£) == 9(x/a,£/a), we get for any 
index m in a bounded range [ 0 ,mo]: 

\\[0p H (9 a )-I]V t qlipm\\L*=O(h M ) t 

uniformly for |t| < T e <, index m G [0, mo] and width a > h £ '^ 3 . 
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According to Prop. 4.9 the ID evolved state tp\ N ^ = e -* t P (JV) / h tp 0 is a linear combination of 
/ + 1 squeezed excited states, plus a remainder 0(h l+1 ~ e ). Hence, taking l = M and using 
the triangle inequality, we get: 

ll[0p,(9„)-/]d A ' l lli==C>(«"). 

We now consider the full state Considering the cutoff 0 as in the statement, we choose 
an auxiliary cutoff 0 (s, t;x, £) = x(t) 0 (x, £), supported near q 0 , such that 0 = 1 near the 
support of 0. We rescale 0 as in the Proposition. If a > h £ / 3 , we observe that in the 
longitudinal variable, we have for H small enough: 

XoihDs)^ 2 ™-^ 3 ) = x((2 im - ip h )h/a){e i{ - 2nn -^ s ) = e *( 2 ™-*«)*. 

As a result, for h small enough we get 

||[Op„(e„) - /idDlf = o(fi M ) 

uniformly for \t\ < T £ > and width a > h £ '/ 3 . 

Let us finally check that the same estimate holds with the cutoff @ a , which is usually not 
of factorized form. From the support proprerties of 0 and 0, we have for any a: 

@q — 1 = (1 — 0a)(9« — 1 ) • 

For a > h £ '/ 3 the symbol calculus in S^ 3 (T x R) shows that this equality translates into 

o Ps (0a) -I = (I- Op ft (0a))(Opn(e Q ) - /) + o(h°°) L ^ L2 . 

□ 

Remark 4.24. In the studies [CR97, Hj99] on the long time evolution of coherent states, 
expansions of the type (4.10) generally break down at the earlier time T\j 3. The breakdown 
is related with the fact that the unstable manifold of the point where the state is centered 
is usually curved in the ambient coordinates ( y , 77 ) used to define the coherent states. Since 
the evolved coherent state wants to expand along this manifold, it wants to curve too, which 
is incompatible with the elliptic shape of squeezed coherent states. This curvature effect be¬ 
comes critical around the time T\/3, leading to the breakdown of the expansion (coefficients 
with large m become dominant). On the opposite, for our normal form Hamiltonian q^ N \ 
the unstable manifold of the fixed point at the origin is the horizontal line, which is not 
curved. The squeezed coherent states can perfectly spread along this line, explaining why 
our expansion does not develop any singularity up to polynomial times t ~ h~ l+e . 

5. L 2 NORMS AND QUASIMODE WIDTHS 

Below we will construct a logarithmic quasimode for the QNF operator p( N \ by averaging 
our evolved coherent states over the time. To begin with, let us compute the energy width 
of the initial Gaussian state i/V 

Lemma 5.1. Consider the initial state given in (4.1), with n G Z possibly In¬ 

dependent, but uniformly bounded when h —> 0 . 

Then if we take the energy 

q° = q°(h(2irn - ^); fr) = E 0 + O(U ), 


(5.2) 
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we obtain the quasimode estimate 

|| (pW _ q°)il) 0 \\ L 2 = \V2 h + o(h 2 ). 

Proof. The separation between the (s, x) variables allows us to replace this norm by a ID 
norm: 

||(p(A) _ g°(^))^ 0 || i2 ( TxR ) = ||Q (jv Vo||l 2 (r) • 

Considering the decomposition = Q^ + Q^P, an explicit computation shows that 

Q[ N) P o = q 1 Op ft ((x£))^ 0 = ihq 1 V^Lp2 , 

Then, Lemma 4.16 shows that each term q a Op^((xf) a )(po in Qn^Vo will produce a linear 
combination of excited states, of norms 0(h a ), with a > 2, so these terms are subdominant 
with respect to the quadratic one. The property q 1 = A + 0{h) achieves the proof. □ 

Hence is a quasimode centered at q°(h) = E 0 + 0(h) and of width ~ Ck. Obviously, 
this is also the case for each evolved state 

5.1. Constructing the logarithmic quasimode for the normal form. In this section 
we will construct a better quasimode for the QNF, using a time averaging procedure. 

Take an arbitrary energy E h = E 0 + 0(h), a time T > 0, a weight function x £ 

C“((—1,1), [0,1]), and its rescaled version xt(£) == x(t/T). Our quasimode is defined 
by 

^xrPn = [ XT(t) e ltEh/h if[ N) dt. (5.3) 

Jr 

It is important to note that this state is not normalized. In order to compute its energy 
width, we will first need to compute its L 2 norm. 

Lemma 5.4. For C > 0 and h e (0,1], we consider a semiclassically large averaging time 
l<T = T h <C\logh\. 

Then the square norm of the state 4 t Xr ,E h takes the form 

IIUt.&II 2 =TS 1 (\, (E t - q°)/h) llxllia (1 + 0(1/T)), 

where 5i(«, •) is a positive function given in (5.12), and q° is the energy (5.2). 

Proof. Like in the previous section, the factorized form of shows that the above norm 
is equal to the L 2 (R)-norm of the ID state 

■Rx.e," f ( 5 . 5 ) 

Jr 

Recall from Prop. 4.9 that the evolved state tp[ N ^ is a combination of squeezed excited 
coherent states: 

21 

Ft dJfqlipo T ^ ^ Cm(t, hJ)TJfql(p2m T Rl | 

m=0 

where the coefficients c m (t,h ) are all 0(h( 1 + \t\ 1 )) and the remainder ||i?/||z ,2 <c(\m l+1 
for C > 0. 
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The square norm of & XT ,E n is expressed by 


1$ 


XT,En I 


GT vA) XT (t') Xr(t) dt dt '. 


(5.6) 


To compute the square norm, we thus need to estimate the overlaps 

(Ttiql(p2m' j Ttql(p2m) = (^2m' j T (t— t')q l ^2rri) ■ 

The first case m — m! — 0 is easy to compute (it is a simple Gaussian integral), and gives 
for any (3 e R [FND03, Eq.40]: 

1 


(ipo,Vp(p Q ) = 


y / cosh(/3) 


(5.7) 


This overlap decreases fast when |/3| 3> 1. We now show that the other terms have a similar 
behaviour. Since Dp = by differentiating with respect to /3 we get 

dpVp = - 1 - Op h {xi)V p = ((a *) 2 ~ a 2 )Vp , 

using the notations of §4 for raising and lowering operators. Since (a*) 2 (po = \Z2(p 2 , we get 

dp{<po,Vp(p 0 ) = V2(tp 0 ,V m ), 

so that 

(ip 0 ,Vp(p 2 ) = 


sinh(/3) 


(2 cosh (a)) 3 / 2 ’ 

Differentiating once more, we obtain similar expressions for (tpo,D a tp^) and (tp 2 , T a (p 2 ). By 
induction, we can obtain in this way explicit expressions for all overlaps ((p 2m ', T a (p 2m ), all 
of which have the form of linear combinations of derivatives of (cosh a)' 1 / 2 . As a result, all 
these overlaps will decay like e~\ a ^ 2 when |a| —» oo. In particular, for any m,m' > 0 there 


exists > 0 such that 


(V^2m' • D a (p 2 m) 


— C m m i , 


uniformly for a G 


(<A), T a (p 0 ) 

(these overlaps are independent of h). 

Expanding the states as in (5.6), the first term (m, m') = (0,0) takes the form 


(5.8) 


/o,o = / {Vo, ©(i-i'H’f’o) e‘ (1 XtV) Xr(t) dt if , where we set 0 h = (E h - q°)/h. 

(5.9) 

From the expression (5.7) we see that the integrand is exponentially localized near the 
diagonal. This motivates us to operate the change of variables t = btk, r — t — t', to get 

r pirdh 


7o,o — 


yj cosh(g 1 r) 


Xr(t + r/2) xr(t ~ r/2) dtdr . 


We are interested in large times T S> 1, so it makes sense to Taylor expand the functions \t 
around the central value t. We use the 1st order Taylor expansion with intermediate value: 


Xt({ ± r/2 ) = Xr(t) ± - X ' T (t ± r ± (t, r)/2) for some r±(f, r) £ (0, r) 


(5.10) 
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The product of the two expansions is split into three terms: 


Xr{t) 2 + ^XT(t)\x' T (t + r+/2) - x '' T (i - r_/2)] - jXt(^ + r+/ 2 )Xr(* - r_/2). 
Let us keep only the first term in the integral I 0 0 , producing 


(5.11) 


Jrd n 


t 1 — 

n — 


dr I xr(t) 2 dt = S^q^ e^T WxWh ■ 


°’° Jr y/cosh(g L r) jr 
T he function <S'i(9 1 ,0r) admits an explicit expression [FND03, Eq.(60)]: 

s,(?1 ' w= ^ |r( i +! 5 r)|2 ' 


2 q 1 


(5.12) 


This function is positive for all values of q 1 ,9n- Given q 1 it takes its maximum at 0% = 0 
with value Si(q 1 ,0) ~ 5.244/g 1 , and decays exponentially when 9 n —> oo. In our situation 
q 1 = A + O(h), and we have 9n = 0(1), so this function is uniformly bounded from below 
and from above by positive constants. 

The second and third terms in (5.11) are both supported in [— T, T] and bounded above 
respectively by y^HxlU 00 Hxlk 00 and ^Hxll! 00 - Injected into the integral Jo,o, these terms 
produce integrals Jq 0 , /g 0 with the following bounds: 


Ko\ < 


T 


l J 0 3 0l < 4T 2 


7 = j dr / = S 2 \\x\\l°°\\x'\\l°o, 

ycosh (q L r) J-t 

[ 7 ==j=^F dr [ ^IIxIU-IIx'IIl- = ^ llx'lli-, 

/ A/cosh(g i r) J-t 1 


where S 2 , S 3 > 0 only depends on q\. When T S> 1, these two terms are subdominant 
compared with Jg 0 . Taking into account that q 1 = A + 0(h), we get 


I 0J0 = TS 1 {\,d H )\\x\\b(l + O(l/T))- 


Let us now consider the parts of involving the corrective terms c m (t, h)T> tq i(p 2 m- 

Such a corrective term may be coupled to the main term T> t / q iip o, or to another corrective 
term h)V t i q np 2 m' , 0 < m! < 21. In both cases, (5.8) shows that the involved scalar 

product T^{t-t') q 1( P 2 m) is bounded above by C m ' m (ip 0 ,'D( t _ t /' )q up 0 ). On the domain of 

integration, all the polynomials are uniformly bounded above by 0[hT l ). These 

terms are thus bounded above by C m . m ' TiT 1 1 0 ${9 h = 0), and are thus much smaller than 
d o,o- 

There remains to treat the contribution to of the remainder Ri. Since all Lp m are 

normalized, a simple computation shows that the parts of ||' 0 ^ 7V ' ) || 2 involving Ri are globally 
bounded above by C h l+1 T l+3 . 

Summing up all contributions, we obtain the stated estimate for the norm of and 

hence the norm of ^> XT ^E h - □ 


Estimating the norm of 4 ’ Xt ,e h was a lengthy, yet necessary step. Using the fact that the 
norm is not very small, it allows us to prove that ^ XT) E h is microlocalized near 70 . Let us 
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define the normalized quasimode 


vh 


def 






XT, Eh 


II ^ Xr,E h || 

Corollary 5.13. If we choose the time T = T e i, the normalized state 4 / Xt ,e h is localized 
in the h £ / 3 neighbourood of 70 , in the sense of Prop, f.22: for any 0 G Cf°{T*(T x R)) with 
0 = 1 in a fixed neighbourhood of 70 , we have the estimate 

ll[Op»(e,.'/i) - /]«W.II^ = C(fi“). (5.14) 


Proof. The proof is a direct consequence of Prop. 4.22. Since 'I' XT ,E h is a linear combination 
of the evolved states we have by the triangular inequality, for any power M > 0: 

ll[Op«(e,.-/i) - < t c .\\ x \\ l ~ o ( h m ) . 

Since the norm of T XTt E h is of order | logfr| 1//2 , we may divide by it on both sides, keeping a 
remainder of the same form. □ 


We have shown in Lemma 5.1 that the states composing ^ XT ,E h are quasimodes for the 
operator P (Ar i(ft) which are centered at the energy q°(h) = E 0 + 0(h). We will take for 
central energy of our quasimode 4 t Xr ,E h the value En used in its definition, which may differ 
from q° by 0{h). To compute the width of this quasimode we need to estimate the norm of 

(pm (n) _ e,)v xtA . 

Proposition 5.15. ForT = T e >, we have the norm estimate 

|| ( pOT _ £»)*„, a || 2 = F S,(A, (E h - q°)/h) ||x'lli= (1 + 0(1/1 logfil)), 

J-e' 

where the function S 1 is defined in (5.12). 

Proof. The proof is a straightforward adaptation of the proof of Lemma 5.4. From the 
definition (5.3) of ^ XT ,E h , we see that 

(pW-E h )* XTiEA = [ X T(t)(P {N) -E h )e-^ p(N) - E ^ 0 dt 

Jr 

= [ XT(t)e itEn/n ihd t fe-^ p(N) - En)/h )i; 0 dt 

Jr 

= -ih [ {d tX T{t))e-^ p(N) - E ^ 0 dt. 

Jr 

The integral in the RHS is similar to the one defining T Xt ,e r ■ To compute its norm, we may 
just apply the Lemma 5.4, up to replacing x by x'■ We get: 

IIT W - B»)*xr,E.H 2 = Ls,( A, (Eh - q °)/fi) HVIlii (1 + 0(1/T)) , 

which is the announced result after specializing to T — XL. □ 

Putting together the results of Lemma 5.4 and Prop. 5.15, we find that the normalized 
quasimodes \k XT) ^ s , centered at the energy E h , have the width 

/(ft) = AMil(l + 0(1/1 log ft|)), 

Te> \\x\\l* 
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which is of the announced order 0{h/\ logfr|). Let us pay attention to the factor. The time 
T e > is the Ehrenfest time (4.21). As in [VS05], we can optimize over the choice of cutoff y to 
Ill'll 2 

minimize the ratio , L . 

M L 2 

Lemma 5.16. 

illf 1ST = 1/2 ’ (5 '17) 

x llxlU 2 

where we take the infimum over all x £ 1), [0,1]). Hence, for any s' > 0, we may 

find a cutoff Xe 1 £ C'^°((—1,1), [0,1]) such that 

Y e '\ L2 — 7r /2(l + s') ■ 
llXe'IU 2 

Proof. By relaxing the condition on the range of y, we may minimize the quadratic form 
q(x) = f(x'(t)) 2 dt over y G C“((—1,1)). This quadratic form can be used to define 
the Laplacian on the interval (—1,1), with Dirichlet boundary conditions: this opera¬ 
tor has the domain Hq((— 1,1)), and admits 1,1)) as a core. Its ground state 

y 0 (t) = cos(7rf/2) reaches the announced infimum. By slightly perturbing y 0 , we 

easily construct a smooth Xe’ with the required property. □ 

Taking the cutoff y £ / in the formula for ^ XTjEn , we get the width 

f(h) = ttA(1 + 2e‘) - ^ ^ (1 + 0( 1/| log^|)). (5.18) 

Provided h is small enough, this width is bouded above by C' 7 , lo l ) /; if we take C 7 = 7rA(l+3e'). 

5.2. A localized quasimode on M: proof of Proposition 1.9. To finish the proof of 
Proposition 1.9, there remains to transform the quasimode ^ xt,e h G S^S 1 x M) for the 
QNF pA) i n to a state fih £ L 2 (M), and check that the latter is a good quasimode for the 
operator P(h). This will be accomplished through the use of the Fourier integral operators 
U E bringing the original Hamiltonian P{h) to the QNF P^ N \ as described in Prop. 3.1. 
Using the notation of this Proposition, we consider the normalized quasimode ^x-rPn with 
the optimized parameters XL, y e / as described in the last subsection, and define the following 
states on L 2 (M), which will be quasimodes for the original operator P(h)\ 

4>h = U N ’L XT ,En ■ 

Here we begin to manipulate our parameters. From Corollary 5.13 we know that ^ X TPn 
up to a remainder 0(h °°), localized in the h £ /2 neighbourood of 7 0 . In this region, the FIO 
U N is essentially unitary, so that 

||t/ ; fi||L 2 (M) = 1 + 0(h°°) . 

Besides, conjugating the localization estimate of the Corollary, we see that V’r is microlocal- 
ized in the h £ /2 neighbourood of 'y(Eo), up to a remainder 0(h°°). 

Now, we want to study the energy width of the quasimode fin with respect to the operator 
P[h). From the conjugation (3.3), we get 

Uj v (P — En) fin = (P^ — ddn)^ XT ,E h + Rn+i^xt,e h ■ 
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The microlocalization of implies the similar localization of (P — E%) ipn, so that 

|| (P - E h ) MlHm) = || U* N (P - E h ) ^|| l2(TxR) + 0(h°°). 

(the FIO U'$t is unitary microlocally near 7 0 , and subunitary away from it). We already 
know the norm of (P^ — E^)^ Xt ^e m given by the optimized width (5.18). There remains 
to check that the norm of Rn+i^ X t,Eh is of smaller order. 

We recall that the symbol of the operator Rn+i has the form tn+i(s, r,x,^,h) = 0((h, x,£) N+1 ) 
when h,x,^ 0. On the other hand, Prop. 5.13 shows that ^ Xt ,e h is localized in a micro¬ 

scopic neighbourhood of y 0 , so we can write 

RN+l^XTtEn = Rn+ 1 Op h {Q h e'/3)^ XT ,E h + . 

The symbol calculus in the class P^p(T*(T x M)) shows that 

R n+ i Op ft (0^v 3 ) = Op fi (r JV+ 1 #0 fie ' /3 ) 

where # denotes the exact Moyal product on T*(T x M). The symbol tn+ i#0 ^'/ 3 is of 
order 0{h e '( Ar+ 1 V 3 ) in the class x M)). The Calderon-Vaillancourt in this class 

then implies that 

llfijr+i o P «(e s ../.)IU^w = o(h^ N + l >' 3 ). 

As a consequence, we obtain 

W(P-E h )^ h \\ LHM) = IKPW - p,)T XTii? j| L2(TxR) +o(n^ N+1 ^ 3 ). 

The first term on the RHS has been estimated in (5.18). To ensure that the remainder is of 
smaller order, we need N, l and e’ to satisfy the condition 

(N + l)£r'/3 > 1, and l> 2. (5.19) 

Remark 5.20. The choice of these parameters should proceed as follows. First we select 
e’ > 0 small, which determines the choice of the time Tb and the near-optimal cutoff \e' 
used to construct our model quasimode. Then, we select N, the order of the normal form, 
large enough to satisfy the condition (5.19): this ensures that the remainder operator Rn+i 
in the QNF, acting on our quasimode ^ Xt ,e h , gives a negligible contribution. The auxiliary 
index / only determines the precision of our expansion for the evolved coherent states, but 
does not impact the definition of the quasimode. 

To end the proof of Proposition 1.9 there remains to check that Vb is microlocalized on 
7 (E 0 ). Using the FIO Un, we may transport the cutoffs Op s (O fie // 3 ) appearing in Prop. 4.22 
and Corollary 5.13, onto the operator 

Un Op,(0,, 7 ,) U’ N = Op,(0) + 0(ft“). 

In each coordinate chart, the symbol O is well-defined up to a remainder in h°°S ( ^(M). 

This symbol is essentially supported in the ^^-neighbourhood of 7 (E 0 ), equal to unity in a 
slightly smaller neighbourhood. The operator Op^(0) is, up to a negligible error, selfadjoint 
(because Un are microlocally unitary near 7). It satisfies 

ll(O P ft(0) — -OV’fill = 0(h °°), 

showing that t/A is microlocalized on j(Eq). 


(5.21) 
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The above microlocalization implies that any semiclassical measure ji sc associated with 
the sequence ( / ipj i ) must be a probability measure supported on 'y(Eo). Because the width 
f{h) = o(h), this measure must be invariant with respect to the Hamiltonian flow As a 
result, we must have p sc = S^e 0 ) which the unique invariant measure supported on the orbit 
'y(Eo). This achieves the proof of Prop. 1.9 

□ 

5.3. Comparison with another quasimode construction. In [CdVP94] Colin de Verdiere 
and Parisse describe 0(h°°) quasimodes and eigenstates of a ID Hamiltonian with a hyper¬ 
bolic fixed point (x, £) = (0, 0). Near the fixed point the Hamiltonian can be put in the QNF 
Q( N] (h). Their strategy to construct the quasimodes is different from the one we presented 
here. Their starting Ansatz are the distributions on <S'(M): 

(pe(x) = f \x\~ 1 ^ 2+ld , for any 9 G M. 

A direct computation shows that (p$ is an eigendistribution of the quadratic operator Q^: 

Q[ N) V t 3 = hq l d<pp. 

From the decomposition (4.3) and Lemma 4.13, <pg is also an eigendistribution of the full 
operator Q^ N \ with eigenvalue hq 1 9 + 0(h 2 ). This distribution appears to be a good starting 
point for a quasimode construction. In [CdVP94] the authors use the fact that the fixed point 
is homoclinic: outside a neighbourhood of the fixed point, the unstable manifold (x-axis) 
bends and becomes the stable manifold (£-axis). This way, they are able to connect together 
the stable and unstable Lagrangian branches of <pp using WKB theory, and produce 0{h °°) 
quasimodes of P(h). 

In our present setting of a “generic” fixed point, we have no information on the continuation 
of the stable and unstable manifolds. The best we can do to the distribution (jpg is to 
microlocalize it inside a compact neighbourhood of (0, 0) G T*M. For this aim we may use a 
cutoff 0 G C£°([—2, 2] 2 ), 0 = 1 in [—1, l] 2 . The state 

<I>e, 0 = Op£(0)w 

is now in L 2 (M), with square norm |||| 2 = 2| log Ti\ + 0(1), each of the four branches of 
the stable and unstable manifolds carrying one fourth of this weight [NV97]. One can show 
that the normalized state Tey = converges to the semiclassical measure <5(o,o)- 

This state is a quasimode of of central energy hq 1 9. Let us compute its width: 

((?<"> - hqtyHes = \Q {N \ Op^e)]^ + Op^(e) (q<"> - ttfe)?, 

= -ifiOp,({g (N) ,e} + 0 (h)) if 

Because 0(x, ^) = 1 near the origin, the symbol {gAO,©} vanishes near the origin and is 
supported in an annulus. As a result, the state [QA9, Op“(0)] is microlocalized along the 
four branches, away from the origin, and has opertor norm 0(h). Finally, we find that $© 
is a quasimode of width 0( i log 1 / 2 )■ This is less sharp than the width 0(K/\ logh|) we have 
obtained by time averaging. 

We now proceed towards understanding the connection between these two quasimode 
constructions. If we replace in the integral (5.5) the smooth time cutoff x by a sharp cutoff 
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H[_i ;1 ], we obtain a quasimode with an energy width 72 ), which resembles the state 

with the identification hq l 9 = En — q°. On the opposite, if we try to represent the 
state & XT ,E h in the form of a truncation Op^(©)<£>e, we obtain a cutoff function of the form 

Q(x, £) ~ x (j outside the disk {( x 2 + £ 2 ) < h}. The corresponding symbol O 

is quite singular when (a: 2 + £ 2 ) > h as it belongs to the “critical” symbol class Spj^ (T*M). 
This singular symbol would be the price to pay if one wants to recover the small energy 
width 0(h/ 1 log/i|) using this construction. 

6. Partially localized quasimodes and small logarithmic widths 
In this section we will prove our main Theorem 1.7. 

We will use the assumption of compactness for M, which ensures that the spectrum of 
P(h) is purely discrete. In Prop. 1.9 we have constructed, for any given sequence ( Eh = 
E 0 + O(h)) n ^ 0 , a family of quasimodes centered at E h and with corresponding 

widths C y | lo g^ , for a constant C 7 = 7 tA (1 + 3e'). In this section we will make the choice 
E h = E 0 for all h. 

For any eigenvalue Ej = Ej(h) e Spec (P(h)), let IT E .ip h be the corresponding spectral 
projection of the quasimode ipn on the 77,-eigenspace. For some C 2 > 0 to be chosen later, 
we consider the spectral interval / = [E 0 — c 2 7 lo ^ ^ , E 0 + c 2 1 1( ^ /t | ] . The Pythagorean theorem 
ensures that 

II {P{h) - E 0 )ijj h \\ 2 = ^ I Ej - E 0 \ 2 \\Il Ej iP h \\ 2 + Y, \ E J - ^run^^ll 2 - (6.1) 

Ej£l c 


Using the quasimode property and the form of the interval I, the second term on the RHS 
above satisfies the following bound, for h G (0, tin)' 

(^y ) 2 * E IU - ^riin^ll 1 > (ji^) V,*ll a ■ (6.2) 

As a result, we get 

l|n^ll 2 < (U) 2 . (6.3) 

To ensure that this upper bound is nontrivial we choose c 2 > C y . As a consequence, the 
projection of iph inside / satisfies 

||n 7 ^|| 2 >l-(^) 2 >0, (6.4) 

in particular the interval / contains at least one eigenvalue Ej of P(h). 

In order to create quasimodes of smaller width than we will project ipn on “short” 
spectral intervals. Precisely, for any choice of e > 0, we take K > 0 large enough such that 


for instance by taking K = [c 2 /e] + 1, 


(6.5) 
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and partition the interval / into K disjoint subintervals (Ik)k=i,...,K of widths ^ and centered 
at the energies e*, = E 0 + , ^7 (—1 + ^j^)- By Pythagoras we have 

l|n r *ll 2 = Xhin 4 « 2 . 

k =1 


For each h , let k(h) be the index (or one of the indices) for which the norm ||ri/ fc '0/ l || is 
maximal. From the lower bound (6.4) this means that 


in, 


k(h) 


The normalized state 


h 


it 


= 


k Q (h) 




( 6 . 6 ) 


(6.7) 


is automatically a quasimode of P(h), centered at En = e^n) and of width e ,^, . 

Let us now study the localization properties of the states -f/y,. The orthogonality of the 
eigenfunctions of P(h) implies the identity ||n/ fco(fi) ^|| 2 = ( 7 / 4 , fph)' 2 - Let us insert the mi¬ 
crolocal cutoff Op^(0) used in (5.21): 

(irni’h ) 2 = A, Op ft(©)Vfi ) 2 + O(ft') 

= (Op «(e)y«,*> 2 + o(ft') 

< II 0p,(e)*|| 2 + o(fi'), 


where we used the fact that Op^(0)* = Op s (0) -\-0(h°°). From these inequalities we deduce 
the lower bound 

II Op,(6)*|| 2 > 4(l - (t^) 2 ) - O(fi'). (6.8) 

This estimate shows that any semiclassical measure /i. sc associated with the sequence of 
quasimodes (ipff) contains a singular component WjS^e 0 ), with a weight 

> ^( 1 - {^r) 2 ) ■ (6-9) 

In other words, the quasimodes (ifn) exhibit a “strong scar” on the orbit 'j(Eq). 

Finally, we may optimize the parameters in the following way: for a given (small) width 
e > 0, we want to keep the weight w 1 as large as possible. By playing a bit with the 
parameters C 2 and K , we end up with the following estimates. 


Lemma 6 .10. For any 0 < e < C 7 , there is a choice of K and c 2 maximizing the RHS of 
(6.9). 

When e <C C 1 , the optimal parameters satisfy the estimates c 2 = y/3C 7 + 0(e), K = K = 
[c 2 /e] + 1, leading to the following lower bound on the weight: 

O^e/C^ 2 ). ( 6 - 11 ) 

Since we can take C y arbitrary close to 7 tA 7 , we get the estimate (1.8) of our Theorem 1.7. 

□ 
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